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We study Kolmogorov-Johnson-MehbAvrami (KJMA) theory of phase conversion in finite vob 
umes. For the conversion time we find the relationship Tcon = Tnu [1 + fd{q)\- Here d is the space 
dimension, Tnu the nucleation time in the volume V, and a scaling function. Its dimensionless 

argument is g = Tcx/tViu, where Tox is an expansion time, defined to be proportional to the diameter 
of the volume divided by expansion speed. We calculate fd{<l) in one, two and three dimensions. 
The often considered limits of phase conversion via either nucleation or spinodal decomposition are 
found to be volume-size dependent concepts, governed by simple power laws for fd{q)- 
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Phase conversions are of importance in physics, chem- 
istry and other fields. Examples are numerous and in- 
clude crystal physics [l|, metallurgy f^ , polymer physics 
@i 13 1 ferroelectric domain switching ^ , magnetization 
and metastability in statistical physics models [l, Q, 
phase transitions in particle physics Q, as well as eco- 
logical landscapes ^] . 

Specific phenomena are nucleation and spinodal de- 
composition Conventionally, for a review see 11 1, 
nucleation is characterized by metastability, while spin- 
odal decomposition is considered to be the mechanism by 
which phase conversion occurs in an unstable system. We 
shall discuss the crossover of these phenomena as func- 
tion of the nucleation time, the expansion speed, and the 
volume. 

In Kolmogorov-Johnson-Mehl-Avrami (KJMA) theory 
13, Hill phase conversion is based on the rate of nucle- 
ation of critical nuclei [l^l and their subsequent expan- 
sion speed due to gain in free energy. Independently, this 
approach was developed a few years later in a concise pa- 
per by Evans [l6| . More recent work derived space-time 
correlation functions 17[ and dealt with screening effects 
[3]. KJMA theory is formulated in an infinite volume. 
But in physics there are no truly infinite volumes. Our 
investigation of finite volumes leads to interesting scaling 
laws. 

There are several scenarios of KJMA theory. Some 
deal also with non-critical nuclei [IJI- We consider here 



the limit in which critical nuclei, compared to the total 
volume are small enough to be considered pointlike. Ex- 
tension of our considerations are possible, but would at 
present distract from the main point. 

In accordance with KJMA theory we make the follow- 
ing assumptions: 

1. Critical nuclei are created with a constant nucle- 
ation rate R — l/{TnuV) at uniformly distributed 



random positions in the volume V. Let us denote 
by the unit volume, and by t"^ the nucleation 
time (average time it takes to create a critical nu- 
cleus in the unit volume). Then the nucleation time 
in the volume V is given by Tm 
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2. Subsequent growth: A nucleus created at time ti 
covers at time t > ti the spherical volume Vi{t) = 
Cd [v {t — ti)]'^, where v is the expansion speed, d 
the space dimension and Cd a dimension dependent 
factor (Ci = 2, Ca = tt, and Cg = Air/S). 

3. The converted volume Vcon{t) is the union of the 
volumes Vi{t) {Vi{t) ~ Q ior t < ti), intersected by 
the total volume V. 

Assumption 1 allows the creation of nuclei in the al- 
ready converted volume Vcon{t)- From assumptions 2 
and 3 it is clear that they do not contribute to phase con- 
version and, therefore, they are not added to the number 
of nuclei in the volume V. Note that KJMA theory of 
phase conversion is kinetic with no details of the respon- 
sible interactions involved. 

The time it takes to transform the bulk system into 
the new phase is the conversion time Tcon- There is 
some arbitrariness in its definition. In essence any con- 
verted volume in the range 0.5 < l^onl'^con)/!^ < 1 , e.g. 
K:on(7'con)/V^ = 0.90, defines a suitable conversion time. 
Only, V^on(''con)/V^ = 1 is not admissible: Tcon will then 
diverge in the infinite volume limit, because due to sta- 
tistical fluctuations some points stay always unconverted 
in an infinite volume. This is well-known in KJMA the- 
ory and even more obvious for systems with fluctuations 
due to interactions. 

For practical reasons we define the conversion time by 
distributing a finite number of trial points uniformly over 
the volume and its boundaries and taking Tcon as the time 
at which all points are first covered by the new phase. 
The number of points is taken to be a constant, inde- 
pendent of the size of the volume. We restrict ourselves 
to cubic volumes of size V — L'^, and choose as trial 



2 



points the sites of a hypercubic lattice that includes the 
2'' corner points of V. Extensions to other geometries are 
straightforward. In particular geometries can be chosen 
to fit actual experimental situations. 

To calculate the average conversion time turns out to 
be easier than one might expect. There are only two 
independent parameters with the dimension of a time, 
T^^ and an expansion time t^x, which we define by 



The functional dependence Tcnn (t°,,, V) is determined 
by a scaling function fdiq) [3 ^ presented below. In- 
stead of Tjjjj we use the nucleation time t^u of the total 
volume V, 

rnu = ^ with A'^ ^ ^ , (2) 
and the scaling function fd{q) is defined by 

— = [l + f4Q)] with (3) 

'^nu '^nu 

The reduction from three variables to one is a mayor 
simplification. That Tcon/''nu depends indeed only on q 
is shown in the following. Natural independent variables 
are r^^, v and L. Using t„u instead of as indepen- 
dent variable is mathematically equivalent. One can then 
define three transformations, which leave q invariant: 
(1) L[ = AiL, v[ = Xiv; (2) = A2L, t^,„, = A2C; 
(3) W3 = X3V, Tg = r^^/As. Combinations of the three 
transformations allow us to create all values t^^, v' and L' 
for which q' — q holds. Now, nuclei at the appropriately 
scaled positions in the volume V' — {L'Y are created 
with the same probabilities as in the initial volume V . 
In case (1) the change in volume (length) is compensated 
by the increase in velocity, so that the conversion time 
stays constant, t[^^^j^ = Tcon- In case (2) the nucleation 
time is scaled, so that the fixed velocity v creates up to 
scaling in size the same patterns as before. Therefore, the 
conversion time scales according to T2 con = ''^2''"con and 
T2,con/'^2,nu = Tcon/Tnu- Similarly, iu case (3) the change 
in velocity is compensated by the change of the nucle- 
ation time, so that the created patterns stay the same 

and 'r;^,con/'^3,nu = Tcon/T-nu bolds. 

In the limit of large expansion speeds (u — > 00, volume 
fixed) we find 

/r""(g)-A,g for g^O, (4) 

where Ad is a dimension and geometry dependent con- 
stant. In this limit, creation of a first critical nucleus 
takes much longer than its subsequent expansion to the 
size of the volume V . Therefore, creation of several criti- 
cal nuclei is unlikely and becomes the time of metasta- 
bility. The conversion time is determined by the farthest 



away corner of the L volume, once the nucleus is cre- 
ated. By integration over the possible positions of the 
nucleus one finds Ai = 0.75, A2 = 1.0704, A3 = 1.315, 
and (for string theorists) = 2.4110. The hmit (g]) 
describes the nucleation scenario of phase conversion. 

At large q, the function fd{q) is up to a multiplica- 
tive constant also analytically determined. Imagine, we 
calculate the conversion time simultaneously on n'^ non- 
interacting systems with identical parameters (nucleation 
time, volume, expansion speed). The conversion time is 
a random variable, which has the same mean value Tcon 
on each system. Let us combine them into one. For 
Tox ^ Tnu the effects due to propagation of phase con- 
version over the boundaries becomes negligible and t'^^^ 
averaged over the combined system is r^on = Tcon- As 
we have q ^ q' = n'^'^-^q and Tnu — ^ = nr'^Tjiu for 
V ~> V ~ n^V, invariance of the conversion time re- 
quires 

/f^°(g) = B,g'^/('^+i) for q^oo. (5) 

This is the limit of spinodal conversion, obtained in vol- 
umes of fixed size for small expansion speeds, w — > 0. 
Many critical nuclei contribute then to the phase con- 
version. For physical parameters r^^, v fixed, and vol- 
ume V — > 00, i.e., A — !■ cx) in Eq. ([2|), the theory al- 
ways describes spinodal decomposition, because q scales 
as g — > A''"'"^ q. 

This is in contradiction to the mean-field approach, 
which leads on infinite volumes to a nucleation region 
with a so called spinodal endpoint 0, [lH. Within the 
more realistic scenario of KJMA theory the spinodal can 
only be an effective concept for finite volumes. In con- 
trast to the comparison with mean-field theory, our re- 
sults are consistent with studies of magnetic field driven 
phase conversion by Rikvold et al. [2^, in which a "dy- 
namical spinodal field" separates the two regimes. 

Let us turn to the general evaluation of fd{q) by Monte 
Carlo (MC) simulations (here not Markov chain MC). 
The implementation of the nucleation process is rela- 
tively straightforward and allows variations of the expan- 
sion speed, and hence g, over many orders of magnitudes. 
This comes, because we have to implement only kinetics 
and no complicating dynamics (for instance, due to in- 
teractions between spins). We use 100 trial points in ID, 
10 X 10 = 100 in 2D, and 5 X 5 X 5 = 125 in 3D. For a 
volume of edge length one this corresponds, in the lattice 
of trial points, to a lattice spacing of 1/99 in ID, 1/9 in 
2D, and 1/4 in 3D. 

The results in ID, 2D and 3D together with the analyt- 
ical g — > and q 00 asymptotic behavior are presented 
in Fig. [1] on a log- log scale. When taking data our step- 
size was a factor two in q. For crosschecks at a few q 
values various combinations of r^^, v and V were used 
that combine to the same q value. 

Performing Gaussian difference tests (e.g., Ref. [211]). 
the first four data are in each case consistent with the 
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FIG. 1: Scaling function fd{q) versus q. 



FIG. 2: Scaling function: End of the small q region (|4|). 



TABLE I: Bd ^ for our trial points from MC simulations. 





ID 


2D 


3D 


min 
Ud 


0.5 


1 


1 


max 
Hd 


4096 


128 


64 


Bd 


2.3285 (60) 


2.0990 (53) 


2.1427 (43) 


nd, Xd(pdf) 


6, 1.24 


7, 1.07 


10, 0.74 



small q approximation (j4|). For q < (7™'"i 9™'" listed in 
table U the data are found to agree with an relative error 
\fd{q) - fT"'^\q)\/ fd{q) < 5% with the analytical small 
q behavior. In the same way they are consistent with the 
large q behavior (O for q > q^^^^, where the Bd values 
listed in table [J arc determined by one-parameter fits to 
Ud data with the largest q values (chi-squared per degree 
of freedom of the fit, Xrf(pdf) = x^/(nd— 1), is also given) . 
In the sense of these approximations we have nucleation 
for q < spinodal decomposition for q > qd^^^, and a 
crossover region in-between. 

In particular, we have in this classification for a 3D 
cubic box Tnu/'T'con < 0.02 for spinodal decomposition and 
Tnu/Tcon > 0.43 for nuclcation. While the coefficients Ad 
and Bd in Eqs. (|4]) and (O depend on the geometry, the 
power laws do not. Therefore, they can be employed to 
characterize the limits universally. 

In Fig. [5] we exhibit the end of the small q region on 
a scale with higher resolution than that of the logarith- 
mic scale of Fig. [TJ Correspondingly the beginning of the 
large q region in 2D and 3D is shown in Fig. [31 In ID 
the asymptotic behavior sets in for considerably larger q 
values. The reason appears to be that there are not yet 
sufficiently many nuclei participating in the phase con- 
version. At the largest q value of Fig. [3] (q = 512) the 
deviation of the ID MC result from its large q asymp- 
totics is still 11%. 

The average number of participating nuclei, Nnucid is 



MC data 3D 
MC data 2D 
_ MC data 1 D 


, — , — , 














Large c 













































































50 100 150 200 250 300 350 400 450 500 

q 

FIG. 3: Scaling function: Beginning of the large q region ^ 
in 2D and 3D. 



smaller than T(.on/'''nu of Eq. ([3]) as nuclei created inside 
an already converted region do not contribute. In Fig. [J] 
we show the ratio R = -/V,iucioi/(''con/''"nu)- For large q it 
approaches R = 0.38 in ID and R = 0.42 in 2D and 3D. 
These numbers are specific to our choice of trial points. 

We continue with illustrations. Changes in physical 
conditions, for instance of the temperature, can influence 
the nuclcation time r^^, the expansion speed v, and the 
volume V. 

Let us assume a constant volume. If the nucleation 
time varies from Tnu fixed expansion veloc- 

ity V, while we stay in the spinodal region q > q^'^^, 
the scaling q ^ q' ^ {'''nu /T^u)q yields for the conver- 
sion time change r^^^ = (r^u/Tnu)^/^'^"'"^^ Tcon- If in the 
same situation the nucleation time Tnu is fixed and the 
expansion speed varies from v —^ v', we find for the 
new conversion time r^^^ = {v/v'Y^'^'^^^^ Tcon- When 
staying in the nucleation region the corresponding equa- 
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FIG. 4: Ratio of contributing nuclei versus q. 



sion (confined — > deconfined) proceeds by spinodal de- 
composition {q > 64)? The answer is t^^ < 5 x 10~^^ 
seconds. This estimate goes up when the expansion speed 
is slower than the speed of light. 

Conclusions. Our equations will need corrections, once 
the critical nuclei can no longer be considered point- 
like, and their size introduces a new dimensional param- 
eter. Further, correlations between nuclei are presently 
neglected and the constant expansion speed of KJMA 
theory may be a too crude approximation for the actual 
dynamical, stochastic expansion process. Nevertheless, 
we think that the scaling laws outlined in this paper are 
at the heart of the distinction between nucleation and 
spinodal regimes of phase transitions. 

We like to thank Per Arne Rikvold and Sachin 
Shanbhag for useful discussions. This work was in part 
supported by the DOE grant DE-FG02-97ER41022. 
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i) Tcon + Ad (Tnu - T^Jq and 
'^con/'^nu = 1 + Ad {v/v') q, respectively. 

Assume a 2D Ising model on a 100 x 100 lattice is 
prepared in its initial state with all spins down. It is then 
simulated by Markov chain Monte Carlo 2l| below the 
critical temperature and with a magnetic field opposite 
to the initial orientation of the spins. For suitable choices 
of temperature and magnetic field the following numbers 
are realistic: (A) Seven nucleation events in one sweep 
with a subsequent expansion speed of 5 lattice spacings 
in 20 sweeps. (B) One nucleation event in 1680 sweeps 
and a subsequent expansion speed of 50 lattice site in 800 
sweeps. A brief calculation puts case (A) with q — 2800 
solidly into the spinodal asymptotics, while with q = 0.95 
case (B) is at the end of nucleation region. Enlarging the 
lattice to 1 000 x 1 000 sites moves case (B) to q — 953 
into the the spinodal region. 

Consider a metastable liquid in a cubic box of size 
(0.1 meter) ^ with a nucleation time of 1 minute in that 
volume and a subsequent explosion-like conversion at the 
speed of 100 km/h. With g = 6 x 10~^ this is deep in 
the nucleation region. This is no longer true if the same 
system is a pool of size (IOmeter)^. Then we are at 
q = 6 000, though the preparation of such a large homo- 
geneous system may in practice be impossible. 

Conversion times of the order of minutes are observed 
in polyethylene crystallization To be definite, let 

Tcon = 180 s. If the nucleation time for the relevant vol- 
ume is < 3.6 s, we would classify the process as spin- 
odal decomposition, and for Tim ^ 77.8 s as nucleation, 
with the crossover region in the range 3.6 s < Tnu < 77.8 s. 

Let us consider the deconfining phase transition Q and 
choose (Ifermi)'^ as the unit volume which defines r^^. 
Suppose the relevant volume at a heavy ion collider is of 
size (lOfermi)'', and that the deconfined phase spreads 
out at the speed of light once a nucleus is created. What 
is the range of nucleation times so that the phase conver- 
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